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Recent experiments with Feshbach-resonant 133 Cs Bose-Einstein condensates have led to unex- 
plained molecule formation: a sudden switch of the magnetic field to its resonance value, followed 
by a finite hold time and another sudden switch to magnetic field values below threshold, converts 
about a third of the initial condensate atoms into molecules. Based on a model of coherent conver- 
sion between an atomic condensate, a molecular condensate, and magnetodissociated noncondensate 
atom pairs of equal and opposite momentum, we find that population trapping is strongly implicated 
as the physical mechanism responsible for molecule formation in switch experiments. 
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Introduction. -Magnetoassociation creates a molecule 
from a pair of colliding atoms when one of the atoms 
spin flips in the presence of a magnetic field tuned near a 
Feshbach resonance [l| . Initial experiments with magne- 
toassociation of Bose-Einstein condensates (BECs) led to 
dramatic losses of condensate atoms in the neighborhood 
of resonance 0, a collapsing condensate with a burst of 
atoms emanating from the remnant condensate 0, in- 
creased losses for decreasing interaction times Q , and co- 
herent oscillations between remnant and burst atoms || • 
Whereas atom-molecule coherence is both necessary and 
sufficient to the explain these observations, the phase 
space densities are insufficient to warrant a claim to 
molecular condensate (MBEC) formation; hence, sub- 
sequent experiments focused on magnetic-field sweeps 
across the Feshbach resonance |fj, which led to the un- 
ambiguous observation of BEC-MBEC conversion Q • 

The backbone of these experiments is coherent evolu- 
tion between an atomic condensate, a molecular conden- 
sate condensate, and magnetodissociated noncondensate 
atoms pairs of equal and opposite momentum 0, . In 
particular, so-called rapid adiabatic passage arises be- 
cause the ground state of the Feshbach system is all 
atoms far above the molecular-dissociation threshold and 
all molecules far below it, so that a slow sweep of the 
magnetic field from one extreme to the other converts 
atoms into diatomic molecules. The colloquial under- 
standing is that a magnetic field sweep directed from be- 
low to above threshold converts an initial atomic conden- 
sate into noncondensate atom pairs of equal and opposite 
momentum, since molecular dissociation is energetically 
favored above threshold, whereas a sweep directed from 
above to below threshold converts the BEC atoms into 
dissociatively stable molecules; i.e., because the dissoci- 
ation channel opens (closes) in a sweep beginning below 
(above) the dissociation threshold. The catch is that the 
so-formed molecules are highly vibrationally excited, so 
that collision-induced relaxation can rather drastically 
limit the lifetime of molecular BECs created in this man- 
ner. 



With this understanding in mind, magnetoassociation 
experiments with 133 Cs condensates .10j have recently 
focused on creating molecules more efficiently For 
sweeps directed from above to below the threshold for 
magnetodissociation, the molecular conversion efficiency 
was found to saturate at about 10% for decreasing sweep 
rates, which is broadly similar to magnetic-field-sweep 
experiments with Rb [a and Na 0. Faced with this in- 
efficiency, Mark et al. [ll| decided to take another tack 
on making molecules: starting from well above thresh- 
old, they abruptly switch the magnetic field value to its 
resonance position, waited a given amount of time, and 
then abruptly switched the magnetic field value to a po- 
sition well below threshold. Somewhat surprisingly, the 
result was a three-fold improvement in the molecular con- 
version efficiency. Moreover, the improved efficiency was 
independent of whether the system was initially above 
or below threshold, as long as the final switch was to 
a magnetic field value well below threshold. Obviously 
these observations defy an understanding in terms of adi- 
abatic conversion of atoms into molecules 0, Q|. The 
purpose of this Letter is demonstrate that molecular for- 
mation in the switch experiments can be explained as co- 
herent population trapping-Rabi/Josephson oscillations 
that damp to a nonzero value-among the coupled atomic 
condensate, molecular condensate, and magnetodissoci- 
ated atoms pairs. 

Atom, Molecule, and Pair Model— Early theories of 
collective magnetoassociation accounted only for the 
condensates, neglecting any and all noncondensate 
modes 01 • However, rogue 0, 0], or unwanted (l4| . 
transitions to noncondensate atom pairs can occur be- 
cause magnetodissociation of a zero-momentum BEC 
molecule need not take the offspring atoms back to the 
zero-momentum atomic condensate, but may just as well 
end up creating two noncondensate atoms with equal- 
and-opposite momenta. Since the collective condensate 
coupling scales like the square root of the magnetic-field 
width of the Feshbach resonance and the magnetodisso- 
ciation rate scales like the width itself, rogue dissociation 
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is expected to play a dominant role in strong magne- 
toassociation. The prime indicator of rogue relevance is 
the density-dependent frequency lj p = Hp 2 / 3 /m [jl Il3|. 
When the collective-enhanced atom-molecule coupling 
satisfies ft > lj p , the role of transitions to nonconden- 
sate atom pairs needs to be carefully considered in atom- 
molecule conversion. 

The purpose of this section is to introduce our minimal- 
yet-realistic model 0, 0] of a joint atom-molecule con- 
densate undergoing rogue dissociation to noncondensate 
atom pairs, in order to determine the mechanism behind 
the unexplained molecule conversion in the Innsbruck ex- 
periments Hence, we assume that N atoms have 
Bose-condensed into the same one-particle state, e.g., a 



plane wave with wave vector k = 0. Magnetoassocia- 
tion then removes two atoms from this state |1), creating 
a molecule in the molecular state |2), which may then 
dissociate into the noncondensate pair state |k, — k). In 
second quantized notation, boson annihilation operators 
for zero-momentum condensate atoms (molecules) and 
noncondensate atoms of momentum k are denoted by 
a (b) and a k . The field-matter interactions that drive 
the atom-molecule transitions are characterized by the 
Rabi frequency K, the binding energy (detuning) of the 
molecule is denoted by 5] also, collisions between conden- 
sate particles are included with a coupling strength \j 
(more on these quantities in a moment). The Hamilto- 
nian for this system reads 
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where e k = h~k 2 /m is the energy of a noncondensate pair, 
and / k describes the wavenumber (energy) dependence 
of the noncondensate coupling. Collisions with noncon- 
densate atoms are neglected for simplicity. 

The corresponding mean-field theory is derived from 
the Heisenberg equations of motion: 
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In Eqs. J2J), we have taken the continuum limit 
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and we also have introduced the rogue and normal den- 
sities, C(e) ee tfi (a k a_ k ) /(27TW 3 / 4 ) and P(e) ee (a k a k ), 

as well as the rogue coupling £ = fl/(2iru^ 4 ). The effects 
of collective enhancement have been included by scal- 
ing the mean-field amplitudes according to x — > x/y/N 
(where x = a,b). The mean-field fractions \a\ 2 and 2\b\ 2 
are then of the order of unity, and the couplings have 
been redefined as Q — ^/Nk and = NXij. 

Finally, the collective-enhanced couplings are given 
explicitly as f2 = [2TThp\aii\A fJ ,AB/m* 1 ] 1 / 2 and A.y = 
2irhpaij I 'mlj = 1,2), where an (ai 2 , a 22 ) denotes 
the zero-field atom-atom (atom-molecule, molecule- 
molecule) scattering length, A M denotes the difference 



between the magnetic moments of the free atom pair and 
the bound molecule, As denotes the magnetic-field width 
of the Feshbach resonance, and m*„ is the reduced mass 
of the ij-th pair; also, the magnetic field dependence of 
the binding energy is given as S = sgn[an]A tl (B — Bq) / h, 
where Bq is the magnetic-field position of resonance. 
Before moving on, we discuss basic numerics. For 
= Aij = P = and £ ^ 0, simple Fourier analysis 
of an initial bound molecular condensate [b(t = 0) = 1] 
gives the below-threshold binding energy as the real and 
negative pole of lj — S — S(o;) + ii] = 0, where r\ = + and 
the molecular self-energy is defined as 
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The simplest energy dependence for the continuum is one 
that obeys the Wigner threshold law up to some abrupt 
cutoff: / 2 (e) = 9(£a/ — e)- The detuning (binding en- 
ergy) then picks up a term S(0) = £ ^/ejw . In principle, 
the cutoff is infinite, and therefore so is the continuum 
shift of the molecular binding energy. To account for this 
divergence, one defines the so-called physical detuning 
5 = 5 — £ 2 ^JTm , which is finite by definition in the limit 
of an infinite cutoff. In practice, any numerical procedure 
employs a finite cutoff, and the finite shift is accounted 
for in exactly the same manner, and the intermediate de- 
tunings in Eqs. @ are taken as physical (renormalized) 
detunings. The number of noncondensate quasicontin- 
uum states and cutoff are chosen to deliver convergence 
and minimize numerical artifacts. 

Results for Sweep and Switch Experiments. -To enable 
a discussion of the 133 Cs sweep and switch experiments, 
we gather explicit parameters from Refs. |lCl 1 1 1| : an = 
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200a , A M = 0.93^o, A B = 2.1 mG, and B = 19.84 G 
[where a (p ) is the Bohr radius (magneton)]. Although 
exact values are unknown, we approximate the effect 
of atom-molecule and molecule-molecule collisions with 
A12 = A22 = An; similarly, we account for collision- 
induced vibrational relaxation by borrowing an imagi- 
nar y s cattering length from Na (see, e.g., Yurovsky et 
al. Q), so that 3[A i2 ] = 3[A 22 ] = px lO^s" 1 . Cru- 
cially, Sl/ojp ~ 2 for p ~ 10 12 cm -3 , so that rogue dissoci- 
ation to noncondensate modes must be treated carefully, 
despite the narrow resonance [lij . 

In order to ensure a uniform magnetic field across 
the atomic sample, the combined optical-dipole and 
magnetic-field-lcvitation trap is turned off, leading to bal- 
listic expansion over the course of the experiments |lll | . 
Nevertheless, we simply assume a peak density consis- 
tent with ~ 10 ms of ballistic expansion, po = 1.1 X 
10 12 cm -3 [ll|, and account for the initial inhomogeneity 
with a local density approximation. This overestimates 
the density and, in turn, the role of vibrational quench- 
ing; on the other hand, the magnetic-field range here 
is truncated, and the sweep rate decreased, which ulti- 
mately shortens the total sweep time and so underesti- 
mates relaxational losses. Of course, these two factors do 
not necessarily cancel each other out, and the ensuing re- 
sults are only qualitative; the aim is not to quantitatively 
model the sweep experiments anyway, but to explain the 
improved efficiency of the switch experiments. 

In the sweep experiments jjjj , the magnetic field is ini- 
tially tuned 0.5 G above the Feshbach resonance thresh- 
old, Bi — Bq + 0.5 G, which is ramped in a linear fash- 
ion, B(t) = Bq — Bt, to a final value 0.5 G below 
threshold, Bf = Bq — 0.5 G. For experimental ramps 
with B < lOG/s the total sweep time is > 100 ms, 
and the conversion efficiency saturates at ca. 10%. For 
numerical ease the B-field range here is truncated to 
B = Bo ± 20As, and we use B = lG/s, for a total 
sweep time of 84 ms. The results are shown in Fig. ^ 
The agreement is surprisingly good, all things consid- 
ered. Also, as the noncondensate population begins to 
build, it is transferred into the molecular BEC; hence, 
creation of molecules instead of dissociated pairs is more 
about coherent conversion than a closing magnetodissoci- 
ation channel (the effect is even more dramatic for slower 
sweeps, as per Ref. 0). 

The relevant frequency scale is f2 0, ^| , so that adia- 
batic following requires 6 /SI < 1. For po and B as above, 
8/SIq « 8; hence, the center of the cloud is still on the 
outside of the adiabatic regime, and the situation is worse 
on the edges of the cloud. Slower sweeps therefore lead 
to more molecules, but allow more time for vibrational 
relaxation, leading to a saturation of the molecular con- 
version efficiency. Nevertheless, it is only necessary to 
satisfy adiabaticity in the immediate vicinity resonance, 
so it should be worth trying to engineer a ramp that is 
"ultrafast" except in the region B ~ Bq ± Ab- In this 
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FIG. 1: Coherent conversion of a Cs atomic condensate 
into a 133 CS2 condensate by sweeping a magnetic field across 
a Feshbach resonance. Here B = lG/s, making for a total 
sweep time of 84 ms, which is, evidently, plenty of time for 
vibrational relaxation to filch ~ 5% of the MBEC population. 



manner, adiabaticity could be satisfied in the vicinity of 
resonance, while minimizing the total sweep time and, 
thus, vibrational relaxation losses, yielding more MBEC. 

Now we turn to the Ref. [ll| switch experiments. Here 
the magnetic field is suddenly switched to B ~ B in 
a time ~ 10 ms, held for a given amount of time, and 
then suddenly switched to a well-below-threshold value 
where the fraction of molecules is determined. Ballistic 
expansion is not really an issue when modeling the switch 
experiments: the switch occurs in a time ~ 10 ms, and 
the hold times are ~ 10 ms themselves, so it is entirely 
reasonable to use p = 1.1 x 10 12 cm -3 throughout; more- 
over, the switch time here is shortened, without chang- 
ing the physics, to hasten the numerics. The results are 
shown in Fig. [21 Whereas not much happens during the 
switch period, during the hold period the system under- 
goes coherent Josephson/Rabi oscillations between the 
atomic condensate, the molecular condensate, and non- 
condensate atom pairs; however, while the couplings re- 
main, the system subsequently freezes into a superposi- 
tion immune to magnetodissociation. The sudden switch 
back then harvests the frozen population. The indepen- 
dence of the whether the initial B-field is above or below 
threshold is seen by comparing Figs. (a,b). More impor- 
tantly, an intriguing feature of the experiments is an ap- 
parent step-wise threshold in the conversion efficiency at 
~ 15 ms, beyond which the molecular efficiency is ~ 30%. 
Due to an inadvertent lag 0], the B-ficld value does 
not reach the vicinity of resonance as expected, so the 
system is still off resonance; along with the fact that 
trapped-population fraction is depends on the detuning, 
as per Figs. (c,d), the lag leads to an apparent thresh- 
old. Similarly, experimental fluctuations in the magnetic 
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FIG. 2: Coherent population trapping in a Feshbach-resonant 
133 Cs Bose-Einstein condensate, following a sudden (0.1ms) 
switch to the vicinity of threshold. In each panel, the to- 
tal fraction is the uppermost (linearly decreasing) solid line, 
whereas the BEC (MBEC, rogue/noncondensate) fraction is 
given by the oscillating solid (dashed, dotted) line. Pop- 
ulation trapping is independent of whether the system is 
initially above (a) or below (b) threshold. However, the 
amount of population trapped does depend on the final mag- 
netic field, as illustrated by (c) Bf — Bo + As/100 and (d) 
Bf — Bo + Afl/10, although it does not really matter whether 
Bf is above or below threshold (not shown). The rogue frac- 
tion in panel (d) is ca. 10~ 3 and is not visible. 



field value ,11] could easily account for the difference be- 
tween the observed (~ 30%) and the calculated (~ 60%) 
molecular fraction for Bf = Bq. 

Conclusions-Given that faster sweeps apparently con- 
vert fewer an fewer atoms into molecules, improving the 
molecular conversion efficiency by speeding up the sweep 
is truly a surprise. The plot thickens with the observation 
of a final 133 Cs2 fraction that is independent of whether 
the initial magnetic field is tuned above or below thresh- 
old (as long as the final switch is to below threshold), 
in contrast with previous experiments 0, Q . Neverthe- 
less, it is known from past work on photodissociation of 
a negative ion |17| that-given a shaped continuum-some 
of the population may stay permanently trapped in the 
bound state. Indeed, the same effect was predicted for 
strong photoassociation [ljj, and since the theory of a 
photoassociation and magnetoassociation resonance are 
formally equivalent, the answer is obvious in hindsight. 
And yet it is still surprising that, given all the possible 
decoherence avenues, coherent population trapping could 



survive. The results herein strongly suggest that it does, 
making coherent population trapping the leading candi- 
date to explain how a fast switch could get the better of 
a slow sweep. 
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